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' Summary. - The strictly isospectral double Darboux method is applied to the quantum 

Taub model in order to generate a one-parameter family of strictly isospectral potentials for this 

■ 

case. The family we build is based on a scattering Wheeler-DeWitt solution first discussed by 
Ryan and collaborators that we slightly modified according to a suggestion due to Dunster. The 



00 



^ \ strictly isospectral Taub potentials possess different (attenuated) scattering states with respect 

> 

00 ' to the original Taub potential. 
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PACS 04.60 - Quantum gravity. 



\0 , 

Q>^ , PACS ll.30.Pb - Supersymmetry. 

bX)' Quantum cosmology and its supersymmetric extension ||l|] are an interesting "laboratory" 



for techniques of current use in nonrelativistic quantum mechanics. One such technique is the 
strictly isospectral double Darboux method (SIDDM), which is a procedure within Witten's 
supersymmetric quantum mechanics |^. 

Previously, we have applied SIDDM to closed, radiation-filled Priedmann-Robertson- Walker 
(FRW) quantum universes, obtaining a one-parameter family of strictly isospectral FRW quan- 
tum potentials and the corresponding wavefunctions 1^]. The Taub minisuperspace model is 
a separable quantum problem and therefore is well suited for SIDDM. It is our purpose in 
this paper to develop the method for the continuous part of the spectrum of the Taub model. 
The motivation for doing this exercise was found in a paper by Sukhatme and collaborators , 
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who obtained bound states in the continuum in quantum mechanics by SIDDM. However, for 
the cosmological Taub case the method does not produce such states as we shall show in the 
following. 

SIDDM is a technique of deleting followed by reinstating an energy level of a one-dimensional 
(ID) potential V{x) by which one can generate a one-parameter family of isospectral potentials 
Viso{x; A), where A is a real, labeling parameter of each member potential in the set. As a matter 
of fact, Khare and Sukhatme |^] were also able to construct multiparameter families of isospectral 
potentials, but we shall be concerned only with the one-parameter case in the following. SIDDM 
is the physical formulation of a mathematical scheme based on the general Riccati solution, which 
has been introduced by Mielnik [0]. The strictly isospectral techniques are well established for 
nodeless bound states, resulting from any pair combination of the well-known Abraham-Moses 
procedure, the Pursey one, and the "supersymmetric" Darboux one 0. Recently, Sukhatme 
and collaborators extended the usage of the double Darboux construction to energy states in 
the continuum region of the spectrum and obtained several families of isospectral potentials with 
bound states in the continuum. In particular, they performed the double Darboux construction 
on the half line free particle wavefunction uq = sin kr, obtaining bound states embedded in the 
continuum for this simple but relevant case. 

Let us take a simple ID Schrodinger equation for a wavefunction u{x) at an arbitrary energy 
level E in the standard form —u (x) + Vb{x)u{x) = Eu{x). Then, the strictly isospectral 
potentials with respect to 14 obtained by means of SIDDM, by which one delets and reinserts 
the level E in the spectrum will be |p 



where T{x) = /^^^ u'^{y)dy and A is the so called isospectral parameter, which is a real quantity. 
The isospectral family of potentials should be understood in the sense that the whole family has 
the same supersymmetric partner "fermionic" potential given by Vf{x) = Vf,(x) — (Inn). 
When u is a nodeless state there are no difficulties in applying the double Darboux technique. 
However, when one uses an eigenfunction with nodes, the corresponding "fermionic" potential 
will have singularities. Nevertheless, the resulting isospectral family of potentials is free of 
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singularities Q and this makes the method viable. Another basic result of the scheme is that 
the u wavefunction is changed into u/[Z{x) + A], i.e., a spatial damping is introduced. 

Let us pass now to the quantum Taub model which has been studied in some detail by Ryan 
and collaborators Q who found that it is separable. Indeed, the Taub Wheeler-DeWitt equation 
is 

52^ 



8a^ gf,.+^"nm = 0, (2) 

where V{f3) = |(e~^^ — 4e~2^). Eq. (2) can be separated in the variables xi = 4a — 8/3 and 
X2 = 4q — 2/3. Thus, one gets two independent ID problems for which the supersymmetric 
procedures are standard practice. The two equations are as follows 

and 

+ -e^2^2 = uj^U2 , (4) 



1 ^2 _ 2 
" dx^o ^ 9 
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where we have already multiplied both sides in Eqs. (3) and (4) by —1 in order to get stan- 
dard Schrodinger equations. The quantity iv is mathematically the separation constant, which 
physically is related to the wavenumber of a positive energy level. 

Martinez and Ryan have considered a wavepacket solution made of wavefunctions ^ having 
the form of a product of modified Bessel functions of imaginary order. We shall slightly modify 
their ^ as follows 

o 3 3 

since, according to Dunster [^, the L function defined as 

L2iu} = ^ . Ahiuj + I-2iuj) , (6) 

zsinh(2Li;7rj 

contrary to the l2iuj function used by Martinez and Ryan, being real on the real axis is a better 
companion for the K function of imaginary order which is also real on the real axis. In order to 
introduce such a change, Dunster invoked the criteria on the choice of standard solutions for a 
homogeneous linear differential equation of the second order due to Miller Q . In the following. 



we shall make the double Darboux construction on the base of a ^' wavefunction of the type 
given in Eq. (5). 

The isospectral potential for the xi variable will be 



1 , 
144 dxl 



Vi{xi; Ai) = — e^'i - 2— ln[Ai + I{xi)] (7) 



and for the X2 one 



where the integrals are 



and 



V2{X2; A2) = le^' - 2-^ ln[A2 + Iix2)] , (8) 



dx 



J(xi) = r Kl{\ey'^)dy (9) 

J-00 

[L2^^{-ey'^) + K2^^{-eyl^)fdy . (10) 
The total Taub isospectral wavefunction has the following form 

.,T , , KU\e^''')L2^.[le^^') + K2^^{leh^') 

^iso{Xl,X2;\l,\2) = Uiso,lUiso,2 = r^, ^ , , n r^. A i X 1 ' 

[Z(xi) + AiJ [Z(x2) + A2j 

We have plotted the Taub isospectral potentials and the corresponding isospectral wavefunction 
for several values of the parameters in Figs. 1, 2, 3, respectively. We confirm the previous findings 
on the strictly isospectral supersymmetric effects Q, except for the fact that the wavefunctions 
are still not normalizable (i.e., they are not bound states in the continuum). 

In conclusion, SIDDM allows the introduction of a one-parameter family of isospectral quan- 
tum Taub potentials having more attenuated states in the continuum region of the spectrum in 
comparison to the original potential. 

The double Darboux method appears to be a quite general and useful method to generate 
new sets of quantum cosmological solutions. This is so because any potential in the Schrodinger 
equation has a classical continuum of positive energy nonnormalizable solutions. Ryan and 
collaborators were among the first to pay attention to the continuum part of the Wheeler- 
DeWitt spectrum. Selecting by means of some preliminary physical arguments that are of 
quantum scattering type one of the continuum solutions one can perform the double Darboux 
construction on that state and generate by this means strictly isospectral families of cosmological 



potentials as well as isospectral cosmological wavefunctions. The quantum Taub model just 
illustrates this nice feature of the double Darboux method. If one pushes further the picture, 
one might say that the incipient universes were nothing else but sets of strictly isospectral states 
of a quantum continuum. The parameter A looks like a decoherence parameter embodying a 
sort of "quantum" cosmological dissipation (or damping) distance. Finally, one can easily apply 
the other procedures of deleting and reinserting energy levels, i.e., combinations of any pairs of 
Abraham-Moses procedure, Pursey's one, and the Darboux one |6|. However, only the double 
Darboux method used here leads to reflection and transmission amplitudes identical to those of 
the original potential. 
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Appendix 

We quote here some important properties of the modified Bessel functions of imaginary order 



Their behavior at x — > 0^ is as follows 



and 
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uj sinh(i^7r) 



vr 



1/2 



[sin(a; ln(x/2) -0) + O(2;2)] 



1/2 



[C0s(lx' ln(x/2) -(/)) + 0(x2)] 



(^1) 



iA2) 



u> sinh(a'vr) 

where the phase (j) = arg[r(l + iiv)]. The amplitudes of oscillation of the two functions become 

unbounded as cj ^ in the neighborhood of the origin. In the paper we worked at fixed, small 

positive CO for computational reasons. 

In the complex plane, as 2; — > oo, the asymptotics are the following 

.1/2 



vr 
2I 



e~^[l + 0(l/z)], |argz| < 3vr/2 - 5 



(^3) 



5 



and 

1 r_nV2 



2i 



e'[l + 0{l/z)], |argz| < 7r/2 - (5 {A4) 



sinh(a;7r) 

where 6 is an arbitrary small positive constant. 

As for the zeros of these functions, it is known that Ki^{x) has an infinite number of simple 
positive zeros in < a; < a; and no zeros in a; < x < oo, whereas Li(^{x) has an infinite number 
of simple positive zeros. Up to x = a; the two sets of zeros are interlaced. 
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Fig. 1 

The original potentials Vi{xi) and ^2(0:2) (left side) and the corresponding isospectral members for 

Ai = A2 = 1 (right side) and w = 1/4. 
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Fig. 

The wave functions ui (the smaller amplitude one) 

Taub universe 
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and U2 (the larger amplitude one) of the original 
for w = 1 /4. 
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Fig. 3 

Isospectral wavefunctions Uiso.i for Ai = 1 (larger amplitude) and Ai = 10 (smaller amplitude); Uiso,2 

for A2 = 1, respectively, and cj = 1/4. 
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Fig. 4 

The original Taub potential in the xi, X2 variables. 
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